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Aims Atrial fibrillation (AF), the commonest cardiac arrhythmia, has been strongly linked with arrhythmogenic sources near the 

pulmonary veins (PVs), but underlying mechanisms are not fully understood. We aim to study the generation and susten- 
ance of wave sources in a model of the PV tissue. 

Methods A previously developed biophysically detailed three-dimensional canine atrial model is applied. Effects of AF-induced 

and results electrical remodelling are introduced based on published experimental data, as changes of ion channel currents 

CcaL. 'ki. 'to. and /Kur). the action potential (AP) and cell-to-cell coupling levels. Pharmacological effects are introduced 
by blocking specific ion channel currents. A combination of electrical heterogeneity (AP tissue gradients of 5-12 ms) 
and anisotropy (conduction velocities of 0.75-1 .25 and 0.21 -0.31 m/s along and transverse to atrial fibres) can results 
in the generation of wave breaks in the PV region. However, a long wavelength (171 mm) prevents the wave breaks 
from developing into re-entry. Electrical remodelling leads to decreases in the AP duration, conduction velocity and wave- 
length (to 49 mm), such that re-entry becomes sustained. Pharmacological effects on the tissue heterogeneity and vulner- 
ability (to wave breaks and re-entry) are quantified to show that drugs that increase the wavelength and stop re-entry (/ki 
and /Kur blockers) can also increase the heterogeneity (AP gradients of 26-27 ms) and the likelihood of wave breaks. 

Conclusion Biophysical modelling reveals large conduction block areas near the PVs, which are due to discontinuous fibre arrange- 

ment enhanced by electrical heterogeneity. Vulnerability to re-entry in such areas can be modulated by pharmacological 

interyentipns. 

Keywords Atrial arrhythmias • Pulmonary veins • Computational modelling • Re-entrant waves • Drug effects 



Introduction 

Atrial fibrillation (AF) self-perpetuates by multifactor electrical re- 
modelling in the atria, which can explain the high prevalence and 
treatment-resistant nature of this commonest sustained arrhythmia. 
The myocardial sleeves of the pulmonary veins (PVs) extending into 
the left atrium (LA) have been recognized as the primary sources of 
high-frequency electrical activity during AF, and PV ablation is rou- 
tinely used to terminate AF clinically.^ However, mechanisms by 



which the PVs can sustain arrhythmogenic activity are unclear, and 
its clinical treatments have significant intrinsic limitations.^ 

Experimental^'^ and computational^''^ studies of the PV sleeves in 
animal models reveal substantial electrical heterogeneities and 
complex arrangements of myocardial fibres resulting in conduction 
discontinuities and unstable re-entrant waves. Similar activation pat- 
terns have been observed at the PV-LA junctions during catheter 
mapping in AF patients,®'' leading to suggestions that high-frequency 
activity in the PV can be sustained by re-entry.'^'^'® Stable re-entrant 
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What's new? 

• Electrophysiologically and anatomically detailed 3D model of 
the pulmonary veins developed to study the effects of electric- 
al tissue remodelling and anti-arrhythmic drugs; 

• Electrical heterogeneity and conduction anisotropy of the 
pulmonary veins both contribute to the generation of 
pro-arrhythmogenic conduction blocks and wave breaks; 

• Cell-to-celi coupling reduction decreases the conduction vel- 
ocity and increases the likelihood of wave breaks, but the 
linked wavelength decrease is not sufficient for sustained 
re-entry; 

• Ion channel remodelling (associated with atrial fibrillation) 
substantially decreases the action potential duration and 
wavelength, and therefore promotes the sustenance of 
re-entry; 

• Potassium channel current blockers can terminate re-entry by 
increasing the action potential duration and wavelength, but 
promote wave breaks by increasing electrical heterogeneity; 

• Amiodarone can both terminate re-entry and reduce the like- 
lihood of wave breaks by substantially decreasing tissue het- 
erogeneity and vulnerability to conduction block. 



patterns have recently been simulated using computational models of 
the human atria:^°the study has demonstrated effects of electrical re- 
modelling on AF self-perpetuation by re-entry, but did not consider 
anisotropy of the PVs. 

Multiple experimental studies have shown that AF-induced elec- 
trical remodelling is characterized by decreases in atrial action poten- 
tial duration (APD),'"" which are associated with underlying 
changes to the expression of several membrane ion channels 
reported both in dog ' and human. Atthe tissue level, remod- 
elling has also been linked with down-regulation and heterogeneous 
expression of connexin proteins that form intercellular gap junctions 
(responsible for the AP conduction), as well as the presence of severe 
fibrosis and fibre disorganization. ^^"^^ The latter factors may con- 
tribute to decreases in the AP conduction velocity (CV)^' and 
increases in conduction anisotropy^ during AF. 

A lack of mechanistic understanding of factors underlying AF onset 
and evolution may explain relatively low success rates in its clinical treat- 
ment, with ~ 30-50% recurrence even after expensive ablation and 
cardioversion procedures.^°''^^ Available anti-arrhythmic drugtherapies 
for AF also have moderate effectiveness and risk of life-threatening 
proarrhythmic complications.^'^ Knowledge of the underlying mechan- 
isms may assist in tailoring the drug therapies to the evolving AF dynam- 
ics, and hence in increasing the efficacy and safety of pharmacological 
treatments for various AF cases. Biophysically detailed models can 
help to dissect such electro-pharmacological mechanisms.'^^"'^^ 

In this study, we apply a detailed electrically heterogeneous and an- 
isotropic model of the canine PV region^ to (i) incorporate effects on 
ion channel and gap junction remodelling seen in chronic AF in 

11 1 2. 1 9 

dog; ■ ■ (ii) investigate the evolution of arrhythmic re-entrant sub- 
strate for AF due to such effects, and (iii) dissect potential pharmaco- 
logical interventions that can terminate the re-entrant activity. 



Methods 

Single cell models 

The Ramirez-Nattel-Courtemanche'^^ model is used to simulate the 
ion channel currents contributing to the AP of a single canine atrial 
myocyte. This model is modified to incorporate detailed the current 
density data specific to the canine LA and PV cells {Figure I and 
Table 1), as in a previous study.'^^ The resultant AP in the PV cell has a 
shorter APD and more positive resting potential than the LA cell 
(figure IB). Atrial fibrillation-induced electrical remodelling is also incor- 
porated into the model based on experimental ion channel current mea- 
surements from dog^'^ (Table 1). This involves a reduction of the 
conductance of the L-type calcium channel current, /cai. by 50%, a reduc- 
tion of the conductance of the transient outward potassium channel 
current, /to, of 55% and an increase in the conductance of the voltage- 
dependant potassium channel current, /ki, of 100%. Such modifications 
are consistent with those observed in human, wherein there are more 
available data.^"'^'*'^^ A previous study has demonstrated the important 
role that remodelling of the ultra-rapid potassium channel current, /^un 
may play in the behaviour of AF in electrophysiologically heterogeneous 
human atrial tissue.^" Because /^ur is atrial-specific, it may be particularly 
relevant to drug studies where it can be targeted without any adverse 
effects on the ventricles. Therefore remodelling of this ion channel is 
also included in the canine cell models. Due to a lack of available canine 
data regarding remodelling of this current, a 50% reduction of /kup 
(Table 1) is incorporated, which is consistent with data from 
human.'^^''^' Ion channel remodelling shortens the APD of both cell 
models (Table 1) and reduces the LA-PV heterogeneity. 

Three-dimensional tissue model 

To simulate electrical activity in the PV at the tissue scale, a previously 
developed three-dimensional (3D) anatomical model of the canine PVs 
and inter-pulmonary region of the posterior LA^ is used (figure JA). 
This model includes fine details of the fibre micro-architecture in this 
region, attained from high-resolution micro-CT images (figure JC and 
D). Details of the reconstruction have been described previously.^'"^" Im- 
portantly for this study, abrupt changes in the fibre arrangements are 
observed at junctions of the PVs and LA, such as the ostium of the left su- 
perior pulmonary vein (LSPV) (Figure ID). Similar abrupt changes have 
been obsen/ed in experimental electro-anatomical studies.'*'^^ 

Action potential conduction is simulated by solving the mono-domain 
equation^^ on a 3D tissue geometry grid, as described previously.^'^"* 
Conduction anisotropy, defined as the ratio between diffusion coeffi- 
cients (characterizing electrotonic cell-to-cell coupling) in the longitudin- 
al and transverse fibre directions, is 16 : 1 in the PVs and 10:1 in the LA. 
These correspond to experimentally validated CV values in the PV and 
LA tissue regions.^ Cell-to-cell coupling reduction is introduced by de- 
creasingthe diffusion coefficients by 50%. This corresponds to the abbre- 
viation in AP conduction velocities associated with AF, experimentally 
observed to be around 10-40% in various stages of remodelling.^^'^^ In 
the 3D model, this resulted in decreases of the maximum longitudinal 
CV from 1.25 to 0.75 m/s and the minimum transverse CV from 0.31 
to 0.21 m/s. Note that atrial tissue remodelling can be underlain by mul- 
tiple pathophysiological factors, such as down-regulation and heteroge- 
neous expression of gap junction proteins, atrial dilatation, and fibrosis. 
However, establishing links with the specific factors was beyond the 
scope of this study. 

Vulnerability to re-entry 

Vulnerability to the initiation and sustenance of re-entrant waves under- 
lying AF is analysed by measuring the time vulnerability windows to 
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Figure I Model of the canine PV region. {A) Three-dimensional segmented geometry, the PV and LA regions are indicated by different colours and 
labelled. (6) Segmentation-specific APs in the LA and PV regions. (C) Fibre orientation in the PV region, atrial fibres are tracked and coloured accord- 
ing to the main local fibre orientation component along the anterior-posterior direction.^ (D) Fibre discontinuities near the LSPV ostium. PV, pul- 
monary veins, LSPV, left superior pulmonary vein; LIPV, left inferior pulmonary vein; RSPV, right superior pulmonary vein; RIPV, right inferior 
pulmonary vein. 



Table I Electrical heterogeneity and remodelling in the PV-LA region 



U 





lc».L(pA/pF) 


/k, (pA/pF) 




/K^pA/pF) 


/K„.(pA/pF) 


APD (ms) 


LA control 


8.0 (8.5) 


5.0 (5.8) 


9.0 (8.4) 


5.0 (5.3) 


10.0 


185 (195) 


PV control 


6.0 (5.9) 


3.5 (3.6) 


4.5 (5.6) 


9.0 (9.8) 


10.0 


173 (175) 


LA remodelled 


4.0 (3.4) 


10.0 (9.3) 


4.0 (3.6) 


5.0 (5.0) 


5.0 


112(125) 


PV remodelled 


3.0 (2.6) 


7.0 (6.1) 


2.0 (2.8) 


9.0 (9.8) 


5.0 


107(130) 



°lon channel current densities and APDs for the PV and 1_A cell models under control and full remodelling conditions. Simulated values are in good agreement with the respective 
experimental values^^ shown in brackets. Note that the current density of /^ur (step current) in control is the same as in the original paper by Ramirez etai. ^The current density of /kuh 
in remodelling conditions is halved as seen in human atrial cells. 



conduction block leading to (i) at least a single wave break (VWcb) and (ii) 
re-entry sustained over several rotations (VWr). An S1-S2 pacing proto- 
col is applied to a site combining substantial electrical heterogeneity (such 
as the LA-PV junction seen in figure J A) and structural anisotropy (such as 
the LSPV region indicated in figure )D), as has been suggested in previous 
studies in canine,^ human^°'^^ and rabbit.^'' Three SI stimuli are applied at 
a cycle length of 280 ms, followed by a short-coupled S2 stimulus. The 
range of S2 stimuli coupling intervals for which an AP can propagate in 



one region (such as the PVs) but is blocl<ed in the other (LA) is considered 
as the VWcB- Within this window, the range of S2 stimuli coupling inter- 
vals for which this conduction block develops into sustained re-entry 
(lasting more than 4 s^°) is considered the VWr. 

Vulnerability windows are assessed in multiple conditions to dissect 
the relative roles of ion channel remodelling, reduced cell-to-cell coup- 
ling, electrical heterogeneity and structural anisotropy in the develop- 
ment of re-entry. The following conditions are considered: control, full 
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Table 2 Drug effects on the APD of the PV and LA cell models in the ion channel remodelling case" 





Control 


Remodelling 


Remodelling + /ki blocl<er 


Remodelling + /Kur blocl<er 


Remodelling + amiodarone 


LA APD (ms) 


185 


112 


189 


138 


92 


PVAPD (ms) 


173 


107 


163 


111 


89 


AAPD (ms) 


12 


5 


26 


27 


3 



°Drug effect on the electrical fieterogeneity between the PV and l_A cells is also quantified as AAPD. 



electrical remodelling (ion channel remodelling and cell-to-cell coupling 
reduction), ion channel remodelling only, cell-to-cell coupling reduction 
only, full remodelling without heterogeneity and full remodelling without 
anisotropy. Heterogeneity is removed by considering all 3D tissue geom- 
etry points as PV cells. As a result, differences in both the APD and resting 
potential between the LA and PV tissue {Figure IB) are eliminated. Anisot- 
ropy is removed by setting the anisotropy ratio to 1 : 1 and adjusting the 
diffusion coefficient such that the overall activation time is comparable to 
control' 

Pharmacological effects 

Models also incorporate the effects of pharmacological ion channel 
blockers, to study their potential effectiveness in the prevention and ter- 
mination of re-entry. Up-regulation of /ki associated with AF-induced re- 
modelling has been shown to be a primary factor in the APD 
abbreviation.^^ Hence, this current provides a natural target for reversing 
the effects of AF-induced electrical remodelling. Due to the atrial-specif ic 
nature of /^ur. this channel also presents an attractive target for a potential 
treatment of AF. Along with /ki and /kup blockade, we also investigate 
effects of amiodarone, an exemplary multi-channel blocker that has 
proved effective in the treatment of AF.^' Effects of /ki and /Kur blockers 
are simulated by reducing their conductances by 50% each. Effects of 
amiodarone are simulated by reducing the conductance of /cai arid /ks 
by 40 and 50%, respectively.^^'"^^ Drug effects on the APD and its hetero- 
geneity between the PV and LA cells are summarized in Table 2. 

To study the relative effectiveness of various pharmacological agents on 
the prevention of re-entry in different AF cases, we quantify their effects 
on the VWcB and VWr in two conditions: (i) full remodelling and 
(ii) cell-to-cell coupling reduction only. In these two cases, the pharmaco- 
logical effects are introduced from the beginning of the respective simula- 
tions. To study the effectiveness of various pharmacological agents in the 
termination of re-entry, their effects are introduced in a separate set of 
simulations after 4 s of sustained re-entrant activity. 




Figure 1 Generation of conduction blocks in the PV region. 
Simulated activation patterns following the LSPV pacing are 
shown for different conditions. {A) Control, (6) full remodelling, 
(C) cell-to-cell coupling reduction only, (D) ion channel remodelling 
only, (E) full remodelling without heterogeneity, (/-) full remodelling 
without anisotropy. All activation times are measured from the time 
of the applied S2 stimulus. Regions coloured in black indicate areas 
in which excitation failed to propagate. The asterisk indicates the lo- 
cation of the applied stimuli. 



Results 

Generation of wave breaks and re-entry 

Pacing in the anisotropic and heterogeneous site near the LSPV 
resulted in various patterns of wave propagation through the PV 
region {Figure 2). These simulations show that all factors considered 
(ion channel remodelling, cell-to-cell coupling reduction, heterogen- 
eity, anisotropy) are important in the initiation and sustenance of 
re-entry. Each of these factors alone can result in conduction block 
and wave breaks, although activation patterns {Figure 2) and VWcb 
(Table 3) are quantitatively different. However, only the combination 
of these factors can result in sustained re-entry (Table 3). In such a 
combination, electrical heterogeneity and anisotropy promote a 



relatively large VWcb. ion channel remodelling results in a decreased 
effective refractory period (ERR, linked to the APD), and coupling re- 
duction results in a decreased CV. The short wavelength (WL = 
ERR X CV) promotes the development of sustained re-entry from 
the initial conduction block and wave breaks. 

In control {Figure 2A), conduction block of a wave spreading from 
the LSPV towards the LA occurs within a large VWcb (18 ms, 
Table 3). However, the relatively fast CV and long wavelength 
(0.78 m/s and 171 mm, respectively) prevent the wave breaks from 
re-entering the LSPV region and generating a complete re-entrant 
circuit. With full remodelling {Figure 26), the same pacing protocol 
can generate sustained re-entry. The combination of remodelled 
electrical heterogeneity (5 ms, Table 2) and anisotropy results in a 
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Table 3 Vulnerability of the PV tissue to re-entry initiation (VWcb) and sustenance (VWr)* 





VWcB (ms) 


VWr (ms) 


WL (mm) 


Control (no remodelling) 


221-239 


- 


171 


Cell-to-cell coupling reduction only 


215-243 


- 


102 


Ion channel remodelling only 


105-115 


- 


81 


Full remodelling (ion channel + coupling) 


103-115 


107-109 


49 


Full remodelling, isotropic tissue 


103-110 




49 


Full remodelling, homogeneous tissue 


103-111 


103 


49 


Full remodelling + /Kur blocker 


111-128 




53 


Full remodelling + /ki blocker 


145-167 




68 


Full remodelling + amiodarone 


97-101 




46 


Cell-to-cell coupling reduction + /Kur blocker 


222-246 




105 


Cell-to-cell coupling reduction + /^i blocker 


335-368 




159 


Cell-to-cell coupling reduction + amiodarone 


167-183 




79 



^The vulnerability windows are quantified under various conditions and under effects of various pharmacological agents. The wavelength is measured as WL = ERP x CV, where ERP 
corresponds to the lower bound of the NAVcb- 



substantial VWcb (7 ms, Table 3), whereas the decreased velocity 
and wavelength (0.47 m/s and 49 mm, respectively) allow for the 
wave breaks to complete and sustain a re-entrant circuit. The 
mechanisms of transition from the initial conduction block to sus- 
tained re-entry (through specific re-entrant conduction pathways) 
are discussed in the next section. 

With cell-to-cell coupling reduction only, re-entry cannot be 
initiated {Figure 2Q. The CV is clearly reduced compared to 
control (see in Figures 2A and C), but the wavelength is not further 
reduced due to he absence of ion channel remodelling (Table 3). 
Therefore, the PVs do not fully recover at the time the wave 
returns to the region of initial conduction block, and re-entry fails. 
With ion channel remodelling only (Figure 2D), sustained re-entry 
cannot be initiated. The short APD (Table 2) — and hence, decreased 
ERP and wavelength (Table 3) — allows for a single re-entrant circuit 
to be completed, but the relatively fast CV due to the intact 
cell-to-cell coupling means the wavelength is still too long to 
permit sustained re-entry. Hence, quick self termination of the re- 
entrant circuit occurs. Note that VWcb is much wider in the 
cell-to-cell coupling reduction case compared to ion channel remod- 
elling (Table 3), and hence the reduced coupling is more favourable 
for the generation of wave breaks (Figure 2Q. However, ion 
channel remodelling is more favourable for re-entry (Figure 26 and 
D). Each component of full electrical remodelling provides more sub- 
strate for re-entry than control. 

Both electrical heterogeneity and anisotropy contribute to the 
generation of conduction block and the propagation pattern promot- 
ing sustained re-entry (Figure 2E and F). Thus, the largest VWcb and 
VWr are observed in the full remodelling condition in the presence of 
both anisotropy and heterogeneity: if either heterogeneity or anisot- 
ropy is removed, VWcb is decreased and VWr largely disappears 
(Table 3). If heterogeneity is removed, the abrupt fibre angle 
change observed at the junction of the PVs and LA by the LSPV 
(Figure 1D) is sufficient to result in a conduction block in this region 
(Figure 2£). Within a very small window (<1 ms. Table 3) this con- 
duction block can develop into sustained re-entry. When anisotropy 



is removed, sustained re-entry cannot be initiated (Figure 2F). Note 
that VWcb is very similar in the anisotropic and homogeneous 
cases (Table 3). Therefore, both the fibre arrangement and electrical 
heterogeneity may play an equally important role in the development 
of wave breaks. 

Transition from wave breaks to re-entry 

Transition from wave breaks (arising from a conduction block near 
the PVs) to stable, sustained re-entry is illustrated in Figure 3. Primar- 
ily, an S2 stimulus applied near the LSPV results in the AP propagating 
within the PVs, but blocked from entering the LA (Figure 3Ai). As 
described in the previous section, the conduction block occurs at 
the LSPV site that combines large electrical heterogeneity and anisot- 
ropy. The wave propagates around the block line to the superior part 
of the LSPV and the inferior part of the LIPV (Figure 3Ai). These direc- 
tions are largely determined by the fibre orientation patterns seen in 
Figure 1C and D. The AP propagates along fibre bundles wrapped 
around the LSPV and LIPV, eventually entering into the LA from 
the posterior and inferior directions (Figure 3Aii). The time needed 
for the waves to propagate around the PVs is sufficient for the LA 
tissue to recover, hence allowing the waves to propagate into this 
region (Figure 3Aiii). The wave then propagates back towards the 
site of initial conduction blocl<, where it re-enters the PVs and 
travels towards the initiation site (Figure 3Aiii and /V). Afterwards, 
the same re-entrant pattern is repeated, with the waves following 
the same conduction pathways around the LSPV and LIPV 
(Figure 36). This transitional re-entrant pattern is observed within 
the first second of simulation (Figure 3G). However, after several 
seconds only a single dominant re-entrant wave is sustained around 
the LIPV (Figure 307), whereas the second wave is suppressed. 

Anti- and pro-arrhythmic drug effects 

Application of each ion channel blocker had a significant effect on the 
wave behaviour. Sustained re-entry initiated in the full remodelling 
condition (Figure 3) was terminated by the application of /^-i (but 
not /Kur) channel blocker or amiodarone after 4 s of simulating the 
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t = 10 ms 



t = 80 ms 



t = 100 ms 



t = 120 ms 




B 



t = 140 ms 



t = 170 ms 



t = 200 ms 



t = 230 ms 



^ ^ 




110 ms 




Figure 3 Generation of sustained re-entry in the PV region. Snapshots of re-entrant w/aves and the resultant activation patterns are shown. (A) and 

(B) Initiation and sustenance of re-entry. Timings indicated are following the S2 stimulus. Arrows indicate direction of wavefront propagation. 

(C) Activation patterns during the transient (i) and sustained (ii) rotations of a stable re-entrant wave. 



stable re-entrant rotation. Sustained re-entry could not be generated 
in the full remodelling condition with the application of either /K^,r 
or /ki channel blockers or amiodarone (Table 3) at the start of 5 s 
simulations. 

However, effects of the pharmacological agents on the generation 
of wave breaks were different. Blocking either /Kur or /ki increased the 
wavelength, but also increased the extent of the VWcb (Table 3). 
Thus, blocking /Kur and /ki increased the VWcb from 12 to 17 and 
22 ms, respectively. Therefore, both these channel blockers dis- 
played anti-arrhythmogenic effects through increasing the wave- 
length and decreasing the likelihood of sustained re-entry, but 
pro-arrhythmogenic effects through increasingthe electrical hetero- 
geneity (Table 2) and likelihood of conduction blocks and wave 
breaks (Table 3). Note that conduction block lines and multiple wave- 
lets have been linked primarily with chronic AF.^''^° 

Amiodarone was the only drug considered without apparent 
pro-arrhythmogenic effects. Apart from preventing the generation 
and sustenance of re-entry (see above), it also substantially decreased 
the VWcb from 12 to 4 ms (Table 3), therefore decreasing the like- 
lihood of wave break generation. This is consistent with the clinically 
observed superiority of amiodarone over other drugs used for AF 
treatment.^ ^'■^^'^^ The underlying electrophysiological mechanism 
can be based on the reduction of APD differences between the PV 
and l_A cells by amiodarone (Table 2). As shown above, such 



electrical heterogeneity is important in the generation of both 
wave breaks (Table 3) and re-entry {Figure 2E). 

Note that the same drugs can result in different wave behaviours in 
different AF cases — thus, drug effects in the cell-to-cell coupling re- 
duction only case result in very different ERPs and VWs (Table 3). 



Discussion 

Simulations of a biophysically detailed 3D model of the canine PV 
region show that a combination of factors is needed to generate re- 
entrant waves: (i) electrical heterogeneity and conduction anisotropy 
can lead to pro-arrhythmogenic conduction block lines and wave 
breaks; (ii) cell-to-cell coupling reduction decreases the CV and 
increases the likelihood of such wave breaks, but does not provide 
a substrate for sustained re-entry; (iii) ion channel remodelling 
decreases the ERP and promotes the sustenance of re-entry; 
(iv) /Kur blockers can prevent the generation of transient re-entry, 
but cannot prevent the generation of wave breaks or terminate sus- 
tained re-entry; (v) /ki blockers can effectively eliminate re-entrant 
substrate by substantially increasing the ERP, but also increase the 
likelihood of pro-arrhythmogenic wave breaks by increasingthe elec- 
trical heterogeneity; and (vi) amiodarone can both eliminate re- 
entrant activity and decrease the likelihood of wave breaks. 
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These results are in agreement with previous experimental and 
computational studies. Blockers of repolarising channel currents 
are known to reduce effectively the re-entrant substrate size by in- 
creasing the APD and wavelength.^'' However, currently available 
anti-arrhythmic drugs are only moderately effective in preventing 
or terminating AF, and may also exert adverse effects.'^ Our study 
may explain such adverse effects by the increases of the VWcb and 
the likelihood of wave breaks in response to blocking 1ki and /Kur- 
Such wave breaks (rather than re-entry) have been linked with 
chronic AF.^''^° This may explain experimental data showing that 
the inhibition of /Kur can actually promote the induction of AF.^^ 
Our simulation results may also provide mechanistic insights into 
the clinically observed superiority of amiodarone over other drugs 
used for AF treatment,^ ^'^^'^^ as amiodarone uniquely reduces the 
likelihood of both re-entry and wave breaks (Table 3). 

Results in this study are consistent with those of a previous com- 
putational study in humans.^" In both studies, ion channel remodelling 
decreased but did not eliminate regional heterogeneity. The pres- 
ence of electrical heterogeneity (Table 2) combined with the short 
wavelengths (Table 3) associated with both ion channel remodelling 
(decreased ERP) and cell-to-cell coupling reduction (decreased CV) 
promoted the development of stable high-frequency re-entrant cir- 
cuits in both studies. Another previous study allowed for the first 
detailed investigation of the role of fibre anisotropy in the develop- 
ment and behaviour of re-entrant circuits in the PVs.^ This study 
built upon the previous resu lts^''° by combining detailed electrical 
heterogeneity and anisotropy in the PV region^ with electrical re- 
modeling,^° as well as the pharmacological effects."^'' The results 
obtained allow us to dissect factors responsible for various aspects 
of wave dynamics underlying AF. 

Interestingly, our new results for the first time clearly demonstrate 
that abrupt changes in fibre angle alone can produce conduction 
blocks, wave breaks and sustained re-entry (Table 3). Moreover, con- 
duction pathways determined by the fibres provide the substrate for 
sustaining re-entry (figure 3). These results are also supported by test 
simulations where the LA (rather than the PV) was fast-paced under 
the full remodelling condition. As a result of such pacing, an area of 
conduction block was observed between the LA and PV regions — 
this area had approximately the same location as the conduction 
block area seen during the PV pacing under same conditions 
[Figure 2), indicating the role of local fibre orientation. The VWcb 
due to the LA pacing (108-116 ms) was smaller than that due to 
the PV pacing (103-115 ms), as in the latter case the anisotropy 
effect was enhanced by the PV-to-LA conduction against the APD 
gradient. Thus, the likelihood of wave breaks was lower during the 
LA pacing than the PV pacing, and no re-entry was observed during 
the LA pacing. More detailed studies are warranted to identify the 
most likely areas of re-entry initiation and sustenance in the 3D atria. 

Limitations 

This study considered the PV region only, as AF is commonly asso- 
ciated with this region.'"'' However, it is important to consider the 
behaviour of re-entry within the PVs in the context of the entire 
atria. This would enable investigation of the evolution of waves 
driven by the PVs at the organ scale, and would also provide a 
much larger substrate for wave propagation and re-entry. This may 
change the effectiveness of the drugs in terminating re-entry in the 



entire atria. For example, /ki block was successful in terminating 
re-entry by increasing the wavelength. Within the small physical sub- 
strate of the PVs, the re-entrant wave did not have sufficient space to 
perpetuate, and hence self terminated. The effect on re-entry per- 
petuation may be different at the larger whole atria scale. A biophy- 
sically detailed model for the entire 3D canine atria is currently being 
developed.'*'^ Note also that the regional PV model enables more 
time efficient study of local tissue properties (such as VWcb)- 
Similar studies using the entire 3D atria models can be very expensive 
computationally. 

It is also important to note a simple nature of the models of drug 
actions. This study did not consider drug binding kinetics or dose de- 
pendencies,''^ as well as potential adverse effects on the other parts on 
the heart.'^'^ Such details were not required to achieve the aims of this 
study. Our simulation results are in good agreement with cell-to-tissue 
level electro-pharmacological data,'^'''^^and may provide useful insights 
into the effects of drugs on wave dynamics underlying AF. Moreover, 
cell-to-cell coupling reduction in this study was uniform, and the CV 
changes were not linked with AF-induced increases in conduction an- 
isotropy and levels of fibrosis.^^"'^ Careful consideration should also 
be given to multiple cell electrophysiology factors (such as relative con- 
tributions of the APD and resting potential to electrical heterogen- 
eity), as well as differences between the single cell and intact tissue 
experimental data used in the models.'^''''''''^ Quantifying all these 
factors and incorporating them into 3D models could provide a big 
challenge for both experimental and computational studies of the 
mechanisms of AF sustenance. 

Conclusion 

The complexity of mechanisms and large variability in pathophysio- 
logical factors underlying AF make it likely that distinctive and often 
several therapeutic options will be needed for individual AF cases. 
In this study, we have applied computational modelling for studying 
the evolution and pharmacological modulation of wave dynamics in 
the PV tissue under conditions of AF-induced remodelling. Simula- 
tion results demonstrated that ion channel remodelling and 
cell-to-cell coupling reduction both reduce the wavelength and 
promote the development of stable re-entrant circuits, while pre- 
serving the substantial electrical heterogeneity and anisotropy to 
create the asymmetries necessary for the initiation of wave breaks 
and re-entry. Elimination of individual factors (such as atrial hetero- 
geneity, anisotropy, ion remodelling, cell-to-cell-coupling reduction) 
can lead to the termination of re-entry, but multiple wave breaks still 
persist. Pharmacological modulation of the electrical heterogeneity 
and vulnerability are also quantified to show that specific drugs that 
eliminate re-entry (an anti-arrhythmic effect) can atthe same time in- 
crease the likelihood of wave breaks (a pro-arrhythmic effect). As 3D 
atrial models incorporate increasingly realistic anatomical and func- 
tional parameters linked to the development of AF, this should 
improve our understanding of electrophysiological mechanisms of 
complex wave behaviours underlying various AF cases, as well as 
their pharmacological prevention. 

Acknowledgements 

The authors acknowledge support from the National Institute for 
Health Research (NIHR) comprehensive Biomedical Research 



Arrhythmogenic wave dynamics in pulmonary veins 



423 



Centre award to Guy's & St Thomas' NHS Foundation Trust in part- 
nership with King's College London and King's College Hospital NHS 
Foundation Trust. 

Conflict of interest: none declared. 

Funding 

This work was funded by a project grant (PG/1 0/69/28524) from the 
British Heart Foundation. 

References 

1. HaTssaguerre M, JaTs P, Shah DC, Takahashi A, Hocini M, Quiniou G etai Spontan- 
eous initiation of atrial fibrillation by ectopic beats originating in the pulmonary 
veins. N Engl j Med 1 998;339:659-66. 

2. Nattel S, Shiroshita-Takeshita A, Brundel BJ, Rivard L. Mechanisms of atrial fibrilla- 
tion: lessons from animal models. Prog Cardiovasc Dis 2005;48:9-28. 

3. Calkins H, Kuck KH, Cappato R, BrugadaJ, Camm AJ, Chen S-A et ai 2012 HRS/ 
EHRA/ECAS expert consensus statement on catheter and surgical ablation of 
atrial fibrillation: recommendations for patient selection, procedural techniques, 
patient management and follow-up, definitions, endpoints, and research trial 
design. Europace 2012;14:528-606. 

4. Hocini M, Ho SY, KawaraT, Linnenbank AC, Potse M, Shah D etai Electrical conduc- 
tion in canine pulmonary veins: electrophysiological and anatomic correlation. Circu- 
lation 2002;105:2442-8. 

5. Arora R, Verheule S, Scott L, Navarrete A, Katari V, Wilson E et al. Arrhythmogenic 
substrate of the pulmonary veins assessed by high-resolution optical mapping. Circu- 
lation 2003;107:1816-21. 

6. Aslanidi OV, Colman MA, Varela M, Zhao J, Smaill BH, HancoxJC etai Heteroge- 
neous and anisotropic integrative model of pulmonary veins: computational study 
of arrhythmogenic substrate for atrial fibrillation, interface Focus 201 3;3:201 20069. 

7. Butters TD, Aslanidi OV, Zhao J, Smaill B, Zhang H. A novel computational sheep 
atria model for the study of atrial fibrillation. Interface Focus 201 3;3:201 20067. 

8. Narayan SM. Krummen DE, Clopton P. Shivkumar K, MIllerJM. Direct or coincidental 
elimination of stable rotors or focal sources may explain successful atrial fibrillation 
ablation: on-treatment analysis of the confirm (conventional ablation for AF v/ith or 
without focal impulse and rotor modulation) trial. jAm Coll Cardiol 2013;62:1 38-47. 

9. Kumagai K, Ogawa M. Noguchi H, YasudaT, Nakashima H, Saku K. Electrophysiolo- 
gic properties of pulmonary veins assessed using a multielectrode basket catheter, 
j Am Coll Cardiol 2004;43:2281 -9. 

10. Colman MA, Aslanidi OV, Kharche S, Boyett MR, Garratt CJ, Hancox JC et ai 
Pro-arrhythmogenic effects of atrial fibrillation induced electrical remodelling- 
insights from 3D virtual human atria. _/ Physiol 201 3;591:4249-72. 

11. Fareh S, Villemaire C, Nattel S. Importance of refractoriness heterogeneity in the 
enhanced vulnerability to atrial fibrillation induction caused by tachycardia-induced 
atrial electrical remodeling. Circulation 1998;98:2202-9. 

12. ChaT-J, Ehrlich JR, Zhang L, Chartier D, Leung TK, Nattel S. Atrial tachycardia re- 
modeling of pulmonary vein cardiomyocytes: comparison with left atrium and po- 
tential relation to arrhythmogenesis. Circulation 2005;111:728-35. 

13. Van Wagoner DR, NerbonneJM. Molecular basis of electrical remodeling in atrial 
fibrillation. _/ Mol Cell Cardiol 2000;32:1 101 - 1 7. 

14. Workman AJ, Kane KA, Rankin AC. The contribution of ionic currents to changes in 
refractoriness of human atrial myocytes associated with chronic atrial fibrillation. 
Cardiovasc Res 2001;52:226-35. 

15. Dobrev D, Ravens U. Remodeling of cardiomyocyte ion channels in human atrial fib- 
rillation. Basic Res Cardiol 2003;98:1 37-48. 

1 6. Ausma j, Wijffels M, Thone F, Wouters L, Allessie M, Borgers M. Structural changes 
of atrial myocardium due to sustained atrial fibrillation in the goat. Circulation 1997; 
96:3157-63. 

17. Van derVelden HMW,JongsmaHJ. Cardiac gap junctions and connexins:their rolein 
atrial fibrillation and potential as therapeutic targets. Cardiovasc Res 2002;54:270— 9. 

1 8. Severs NJ, Coppen SR, Dupont E, Yeh H-1, Ko Y-S, Matsushita T. Gap junction altera- 
tions in human cardiac disease. Cardiovasc Res 2004;62:368-77. 

19. Gaspo R, Bosch RF, Talajic M, Nattel S. Functional mechanisms underlying 
tachycardia-induced sustained atrial fibrillation in a chronic dog model. Circulation 
1997;96:4027-35. 

20. Camm AJ, Kirchhof P, Lip GYH, Schotten U, Savelieva 1, Ernst S et ai Guidelines for 
the management of atrial fibrillation: the task force for the management of atrial fib- 
rillation of the European Society of Cardiology (ESC). Europace 2010;12:1 360-420. 

21 . Weerasooriya R, Khairy P, Litalien J, Made L, Hocini M, Sacher F et ai Catheter ab- 
lation for atrial fibrillation: are results maintained at 5 years of follow-up? j Am Coll 
Cardiol 2011;57:160-6. 



22. Dobrev D, Nattel S. New antiarrhythmic drugs for treatment of atrial fibrillation. 
Lancet 2010;375:1212-23. 

23. Noble D. Computational models of the heartand their use in assessing the actions of 
drugs.; Pharmacol Sci 2008;107:107-17. 

24. Workman AJ, Smith GL, Rankin AC. Mechanisms of termination and prevention of 
atrial fibrillation by drug therapy. Pharmacol Ther 201 1;1 31:221 -41. 

25. Zemzemi N, Bernabeu MO, SaizJ, Cooper J, Pathmanathan P, Mirams GR eta/. Com- 
putational assessment of drug-tnduced effects on the electrocardiogram: from ion 
channel to body surface potentials. Br j Pharmacol 201 3;1 68:71 8—33. 

26. Ramirez Rj, Nattel S, Courtemanche M. Mathematical analysis of canine atrial action 
potentials: rate, regional factors, and electrical remodeling. Am j Physiol l-leart Circ 
Physiol 2000;279:H1767-85. 

27. Aslanidi OV. Butters TD, Ren CX, RyecroftG, Zhang H. Electrophysiological models 
for the heterogeneous canine atria: computational platform for studying rapid atrial 
arrhythmias. ConfProc IEEE EngMed Biol Soc 201 1 ;201 1:1 693-6. 

28. Caballero R, de la Fuente MG, Gomez R, Barana A, Amoros I, Dolz-Gaiton P et ai In 
humans, chronic atrial fibrillation decreases the transient outward current and ultra- 
rapid component of the delayed rectifier current differentially on each atria and 
increases the slow component of the delayed rectifier current in both, j Am Coll 
Cardiol 2010;55:2346-54. 

29. Christ T,BoknikP,W6hKS,WettwerE, Graf EM, Bosch RFeto(.L-typeCa^+ current 
down regulation in chronic human atrial fibrillation is associated with increased activ- 
ity of protein phosphatases. Circulation 2004;1 10:2651 -7. 

30. Aslanidi OV, Nikolaidou T, Zhao J, Smaill BH, Gilbert SH, Holden AV et ai Applica- 
tion of micro-computed tomography with iodine staining to cardiac imaging, seg- 
mentation, and computational model development. /£££ Trans Med imaging 2013; 
32:8-17. 

31. Ho SY, Cabrera JA, Sanchez-Quintana D. Left atrial anatomy revisited. CircArrhythm 
Elearophysiol 2012;5:220-8. 

32. Clayton RH, Bernus O, Cherry EM. Dierckx H, Fenton FH, Mirabella L etai Models of 
cardiac tissue electrophysiology: progress, challenges and open questions. Prog 
Biophys Mol Biol 201 1 ;104:22-48. 

33. Aslanidi OV, Colman MA, Stott J. Dobrzynski H, Boyett MR, Holden AV et ai 3D 
virtual human atria: a computational platform for studying clinical atrial fibrillation. 
Prog Biophys Mol Biol 201 1 ; 1 07:1 56-68. 

34. Aslanidi OV, Boyett MR, Dobrzynski H, Li J, Zhang H. Mechanisms of transition from 
normal to reentrant electrical activity in a model of rabbit atrial tissue: interaction of 
tissue heterogeneity and anisotropy. Biophys j 2009;96:798-817. 

35. Zhang H, Garratt CJ,ZhuJ. Holden AV. Roleof up-regulation of/Ki in action potential 
shortening associated with atrial fibrillation in humans. Cardiovasc Res 2005;66: 
493-502. 

36. Shinagawa K, Shiroshita-Takeshita A, Schram G, Nattel S. Effects of antiarrhythmic 
drugs on fibrillation in the remodeled atrium insights into the mechanism of the su- 
perior efficacy of amiodarone. Circulation 2003;107:1440-6. 

37. Benson AP. Aslanidi OV, Zhang H, Holden AV. The canine virtual ventricular wall: a 
platform for dissecting pharmacological effects on propagation and arrhythmogen- 
esis. Prog Biophys Mol Biol 2008;96:187-208. 

38. Holden AV, Aslanidi OV, Benson AP, Clayton RH, Halley G, Li P etai The virtual ven- 
tricular wall: a tool for exploring cardiac propagation and arrhythmogenesis. _/ Biol 
Ph/s 2006;32:355-68. 

39. Allessie MA, de Groot NMS, Houben RPM. Schotten U, Boersma E, Smeets JL et ai 
Electropathological substrate of long-standing persistent atrial fibrillation in patients 
with structural heart disease: longitudinal dissociation. CircArrhythm Electrophysiol 
2010;3:606-15. 

40. Verheule S.Tuyls E, van Hunnik A, Kuiper M, Schotten U, Allessie M. Fibrillatory con- 
duction in the atrial free walls of goats in persistent and permanent atrial fibrillation. 
CircArrhythm Electrophysiol 2010;3:590-9. 

41. Burashnikov A. Antzelevitch C. Can inhibition of IKur promote atrial fibrillation? 
Heart Rhythm 2008;5:1 304-9. 

42. Varela M, Zhao J, Aslanidi OV. Determination of atrial myofibre orientation using 
structure tensor analysis for biophysical modelling. Lect Note Comput Sci 2013; 
7954:425-32. 

43. Davies MR, Mistry HB, Hussein L, Pollard CE, Valentin J-P, Swinton J eta/. An insilico 
canine cardiac midmyocardial action potential duration model as atool forearly drug 
safety assessment. Am j Physiol hieart Circ Physiol 201 2;302:H 1466-80. 

44. Miyauchi Y, Hayashi H, Miyauchi M, Okuyama Y, Mandel WJ, Chen PS et ai Hetero- 
geneous pulmonary vein myocardial cell repolarization implications for reentry and 
triggered activity, i-leart Rhythm 2005;2:1339-45. 

45. Ehrlich JR, Cha TJ, Zhang L, Chartiier D, Melnyk P, Hohnloser SH et ai Cellular elec- 
trophysiology of canine pulmonary vein cardiomyocytes: action potential and ionic 
current properiiies.jPhys/o/ 2003;551:801-13. 



